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Abstract 

Simulation of a system consisted of free particle bouncing on a 
vertically vibrated based is performed. Two different states, which are 
steady and unsteady energy transfer state are found. The vibrating 
based is hold at constant vibration frequency / = 0.1 as the vibration 
amplitude A varied. Sinusoidal form is used. Granular temperature 
T g as function of based velocity and coefficient of restitution is used 
but shown no role in determining energy transfer state of the system. 
Peak of free particle trajectory x m around value 20 seperate region of 
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the two states. 
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1 Introduction 

System of a bouncing particle on a vertically vibrated based is not so simple 
as it sounds and it can be seen as part of more complex system of granular 
materials. Granular temperature of the system is scaled to based velocity 
and coefficient of restitution pp. In randomly vibrating based it exibits an 
inelastic collapse [2] . Study of stability of existence of periodic modes is also 
conducted [3] , which triggered further behavious characterization numerically 
[1]. Completely inelastic particle is also subject of similar system [5J. A 
steady bounce mode and unsteady one are observed in simulation, which is 
explained by term of granular temperature. 

2 Simulation 

A system consists of two granular particles is used in this work. A particle 
named wall particle will oscilate with certain frequency / and amplitude A. 
Other particle will be fallen from certain height h and then collide the wall 
particle and bounce back. The collision and bouncing happen several times 
and then observed in simulation for a time duration. Both particles have the 
same mass m and diameter D. 

A molecular dynamics method (MD) implementing Gear predictor-corrector 
algorithm [6] of 5th order is used in the simulation. Between wall particle 
and free particle a linear spring- dashpot model is used [7J. 

Parameters used in the simulation is as followed: xq — 0, X\ — 4, m — 1, 
d = l } g = 0.5, / = 0.1, ^(0) = 0, k = 100, 7 = 0.1, At = 10" 2 , t { = 0, and 
tf = 1000. The amplitude A is then varied between 0.03125 and 1.8, which 
later represented in granular temperature T g instead of amplitude A. 

3 Results and discussion 

As written in [7] , there is a relation 
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Figure 1: A system consisted of wall particle (left) and free particle (right). 
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where is this simulation it it used that k = 100, 7 = 0.1, and m = 1, so 
it will be obtained that 
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which in this case it is used that 0.03125 < A < 1.8. Left side of Equation 
is also known as normalized acceleration IS]. 
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It is used that g = 0.5, it can be calculated that 0.0247 < T < 1.4212. 
Then right side of Equation gives 0.170, which tells us that the conducted 
simulation does not fit the proposed requirement. Even the requirement does 
not meet it can still be calculated the granular temperature PQ. 
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Using Equation(TfJ) a plot of x m againts T g can be obtained, which is in 
our case 



T g = 213.6A 
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Figure 2: Example of trajectory of bouncing particle at / = 0.1 for differ- 
ent amplitude and state: A = 1.43, unsteady (top) and A = 1.50, steady 
(bottom). 

Based on Figure [3] it can be said that granular temperature T g does not 
play significance role in determining whether the free particle will have a 
steady or unsteady energy transfer state. These states can be differenced 
using maximum peak of free particle trajectory x m , which is around 20 as 
shown in Figure [3j 

4 Conclusions 

Steady and unsteady energy transfer states are defined and can be loosly 
determined using maximum peak of free particle trajectory. The granular 



4 



100 
80 
60 
40 
20 




Q 



0:: / 



O 



unsteady 
steady 

C©C O •<•') o 



.o-e-G- &o°-q- 



50 100 150 200 250 300 350 400 

Figure 3: Region of steady and unsteady state (seperated by dashed line) as 
function of maximum bouncing height x m and granular temperatur T g . 

temperature shows no significance in driving the state in which the system 
stays. 
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